Itinerant in-plane magnetic fluctuations and many-body correlations in Na x .CoC>2 
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Based on the ab-initio band structure for Na^CoCh we derive the single-electron energies and 
the effective tight-binding description for the i2 9 bands using projection procedure. Due to the 
presence of the next-nearest-neighbor hoppings a local minimum in the electronic dispersion close to 
the r point of the first Brillouin zone forms. Correspondingly, in addition to a large Fermi surface 
an electron pocket close to the T point emerges at high doping concentrations. The latter yields 
the new scattering channel resulting in a peak structure of the itinerant magnetic susceptibility at 
small momenta. This indicates dominant itinerant in-plane ferromagnetic fluctuations above certain 
critical concentration x m , in agreement with neutron scattering data. Below x m the magnetic 
susceptibility shows a tendency towards the antiferromagnetic fluctuations. We further analyze 
the many-body effects on the electronic and magnetic excitations using various approximations 
applicable for different U/t ratio. 

PACS numbers: 74.70.-b; 71.10.-w; 75.40.Cx; 31.15.Ar 



I. INTRODUCTION 



The recent discovery of the superconductivity in 
hydrated lamellar cobaltate Na^CoO^ • yEbGi raised 
tremendous interest in the nature and symmetry of the 
superconductive pairing in these materials. The phase 
diagram of this compound, with varying electron dop- 
ing concentration x and water intercalation y, is rich and 
complicated; in addition to superconductivity, it exhibits 
magnetic and charge orders, and some other structural 
transitions^^ 5 -. The parent compound, Na^CoC^, is 
a quasi-two-dimensional system with Co in C0O2 layers 
forming a triangular lattice where the Co-Co in-plane 
distance is two times smaller than the inter-plane one. 
Na ions reside between the C0O2 layers and donate ad- 
ditional x electrons to the layer, lowering the Co valence 
from Co 4+ (3d 5 configuration) to Co 3+ (3d 6 configura- 
tion) upon changing x from (C0O2) to 1 (NaCoC^). 
The hole in the d-orbital occupies one of the ti g levels, 
which are lower than e g levels by about 2 eV— . The de- 
generacy of the t2 g levels is partially lifted by the trigonal 
crystal field distortion which splits the former into the 
higher lying a\ g singlet and the lower two e' g states. 

First principles LDA (local density approximation) and 
LDA+U band structure calculations predict Na^Co02 
to have a large Fermi surface (FS) centered around the 
r = (0,0,0) point with mainly a\ g character and six 
hole pockets near the K = (0,4^,0) points of the hexag- 
onal Brillouin zone of mostly e' g character for a wide 
range of x^. At the same time, recent surface sensitive 
Angle-Resolved Photo-Emission Spectroscopy (ARPES) 
experiments^ ' 10 : 11 reveals a doping dependent evolution 
of the Fermi surface, which shows no sign of hole pockets 
for a wide range of Na concentrations, i.e. (0.3 < x < 
0.8). Instead, the Fermi surface is observed to be cen- 



tered around the T point and to have mostly a\ g char- 
acter. Furthermore, a dispersion of the valence band is 
measured which is only half of that calculated within the 
LDA. This indicates the importance of the electronic cor- 
relations in Na^CoC^. 

Shubnikov-de Haas effect measurements revealed two 
well-defined frequencies in Nao.3Co02, suggesting either 
the existence of Na superstructures or the presence of the 
e' g pockets 1 ^. Last possibility was found to be incompat- 
ible with existing specific heat data. Also, within the 
LDA scheme the Na disorder was shown to destroy the 
small e' g pockets in Na 2 /3Co02 because of their tendency 
towards the localization^. 

The hole pockets are absent in the LSDA+U 
(Local Spin Density Approximation + Hubbard U) 
calculations^. However, in this approach the insulating 
gap is formed by a splitting of the local single-electron 
states due to spin polarization, resulting in a spin polar- 
ized Fermi surface with an area twice as larger as ARPES 
observes. 

The dynamical character of the strong electron corre- 
lations has been taken into account within Dynamical 
Mean Field Theory (DMFT) calculations^ 5 - and, surpris- 
ingly has led to an enhancement of the area of the small 
Fermi surface pockets, in contrast to the experimental 
observations. At the same time, the use of the strong- 
coupling Gutzwillcr approximation within the multior- 
bital Hubbard model with fitting parameters^ yields an 
absence of the hole pockets at the Fermi surface. Ac- 
cording to these findings, the bands crossing the Fermi 
surface have a\ g character. 

Concerning the magnetic properties, LSDA predicts 
Na x Co02 to have a weak intra-plane itinerant ferro- 
magnetic (FM) state for nearly all Na concentrations, 
0.3 < x < 0.7™. On contrary, neutron scattering 



finds A-type antiferromagnetic order at T m w 22K 
with an inter(intra)-plane exchange constant J c ( a b) = 
12(— 6)mcV and with ferromagnetic ordering within Co- 
layer only for 0.75 < x < 0.9±&±&2°. 

In this paper we derive an effective low-energy model 
describing the bands crossing the Fermi level on the basis 
of the LDA band structure calculations. Due to the FS 
topology, inferred from LDA band structure, the mag- 
netic susceptibility xo (q, ui = 0) reveals two different 
regimes. For x < 0.56 the susceptibility shows pro- 
nounced peaks at the antiferromagnetic (AFM) wave vec- 
tor Q AFM = {(f ,(x' )} resulting in a ten- 
dency towards in-plane 120° AFM order. For for x > 
0.56 the susceptibility is peaked at small momenta near 
Qfm = (0, 0). This clearly demonstrates the tendency of 
the system towards an itinerant in-plane FM state. We 
find that the formation of the electron pocket around the 
r point is crucial for the in-plane FM ordering at high 
doping concentrations. We further analyze the role of 
the many-body effects calculated within the Fluctuation- 
Exchange (FLEX), Gutzwillcr, and Hubbard-I approxi- 
mations. 

The paper is organized as follows. In Section [II] the 
LDA band structure and tight-binding model parame- 
ters derivation are described. The doping dependent 
evolution of the magnetic susceptibility within the tight- 
binding model is presented in Section IIII1 The role of 
strong electron correlations is analyzed in Section IIVI 
The last Section summarizes our study. 



II. TIGHT-BINDING MODEL 

The band structure of Nao.6iCo02 was obtained 
within the LDA^i in the framework of TB-LMTO- 
ASA (Tight Binding approach to the Linear Muffin-Tin- 
Orbitals using Atomic Sphere Approximation)^ com- 
putation scheme. This compound crystallizes at 12K 
in the hexagonal structure (P^s/mmc symmetry group) 
with a=2.83176Aand c=10.84312A3. A displacement 
of Na atoms from its ideal sites 2d (1/3,2/3,3/4) on 
about 0.2Aare observed in defected cobaltates for both 
room and low temperatures. This is probably due to 
the repulsion of the randomly distributed Na atoms, 
locally violating hexagonal symmetry^. In this study 
Na atoms were shifted back to the high symmetry 2d 
sites. Oxygen was situated in the high-symmetry 2d- 
position (1/2,2/3,3/4). The obtained Co-0 distance is 
1.9066Awhich agrees well with experimentally observed 
one 1.9072(4)A22. This unit cell was used for all doping 
concentrations. The effect of the doping was taken into 
account within the virtual crystal approximation where 
each Co site has six nearest neighbor virtual atoms with 
fractional number of valence electrons x and a core charge 
10 + x instead of randomly located Na. Note, that all 
core states of the virtual atom are left unchanged and 
corresponds to Na ones. We have chosen Co 4s,4p, 3d 
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FIG. 1: (color online) Calculated near- Fermi level LDA 
band structure and partial density of states (PDOS) for 
Nao.33Co02- The contribution of Co-oi 9 states is denoted by 
the vertical broadening (in red) of the bands with thickness 
proportional to the weight of the contribution. The crosses 
indicate the dispersion of the bands obtained by projection on 
the tig orbitals. The horizontal line at zero energy denotes 



the Fermi level. 
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FIG. 2: (color online) (a) Schematic crystal structure of the 
Co- layer in Naa;Co02 with hopping notations within the first 
three coordination spheres (C.S.). (b) LDA-calculated Fermi 
surface with cylindrical part (violet) having mostly a\ g char- 
acter and six hole pockets (red) having mostly e' g character. 
k x and k y coordinates of the symmetry points are given in 
units of 2-rv/a with a being the in-plane lattice constant. 



states, 2s, 2p, 3d states of O and 3s, 3p, 3d states of Na as 
the valence states for the TB-LMTO-ASA computation 
scheme. The radii of atomic spheres are 1.99 a.u. for Co, 
1.61 a.u. for oxygen, and 2.68 a.u. for Na. Two classes of 
empty spheres (pseudo-atoms without core states) were 
added in order to fill the unit cell volume. 

In order to find an appropriate basis the occupation 
matrix was diagonalized and its cigenfunctions were used 
as the new local orbitals. This procedure takes into ac- 
count the real distortion of the crystal structure. The 
new orbitals are not pure trigonal a\ g and e' g orbitals 
but we still use the former notations for the sake of sim- 
plicity. 288 fc-points in the whole Brillouin zone were 
used for the band structure calculations (12x12x2 mesh 
for k X7 k y , and k z , correspondingly). 

The bands crossing the Fermi level are shown in 
(Fig. [T|). One sees that they have mostly ai g character, 
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consistent with previous LDA findings^. Note, the small 
FS pockets near the K point with e^-symmetry present at 
x = 0.33 [see Fig. |2fb)] disappear for higher doping con- 
centrations because of the corresponding bands sink be- 
low the Fermi level. The difference in the dispersion along 
K — M and L — H directions is due to non-negligible in- 
teraction between C0O2 planes. A small gap between Co- 
3d and 0-2p states at about -1.25 eV present for x=0.61 
disappears for x=0.33 due to the shift of the <i-band to 
lower energy upon decreasing number of electrons. 

In the following we restrict ourselves to the model with 
the in-plane hoppings inside C0O2 layer to describe the 
doping dependence of the itinerant in-plane magnetic 
order. Hence we neglect bonding-antibonding (bilayer) 
splitting present in the LDA-bands. This assumption 
seems to be justified since the largest interlayer hopping 
matrix clement is an order of magnitude smaller than the 
intra-planc one (0.012 cV vs. 0.12 eV). 

To construct the effective Hamiltonian and to derive 
the effective Co-Co hopping integrals t°f g for the t 2g - 
manifold we apply the projection procedur e) 24 ' 25 . Here, 
(a/3) denotes a pair of orbitals, a\ g , e' gl or e' g2 . The 
indices / and g correspond to the Co-sites on the tri- 
angular lattice. The obtained hoppings and the single- 
electron energies are given in Table U for three difference 
doping concentrations. A comparison between the bands 
obtained using projection procedure and the LDA bands 
is shown in Fig. [T] confirming the Co-t 2 g nature of the 
near-Fermi level bands&2£. For simplicity we have enu- 
merated site pairs with n = 0,1,2,..., t°f g -> (see 
Fig. [2](a) and the correspondence between in-plane vec- 
tors and index n in Table [I]). Due to the C3 symmetry of 
the lattice, the following equalities apply: {t^ 13 ] = 
\tf\ = \tf\, \tf\ = \tf\. In addition tf = tf for 
dig — ► dig hoppings, which, however, does not hold for 
e g i 2 orbitals. Since the hybridization between the a\ g 
and the e' g bands is not small, a simplified description 
of the bands crossing the Fermi level in terms of the a\ g 
band only (neglecting e' g band and the corresponding hy- 
bridizations, see for example Refill) may lead to an in- 
correct result due to a higher symmetry of the ai 9 -band. 

In summary, the free electron Hamiltonian for C0O2- 
plane in a hole representation is given by: 

Ho = -J2 ( £ ° »W "EE *f d L<A/^ (1) 

k,a,<7 k,<r a,/3 

where d^ aa (d\^ a<y ) is the annihilation (creation) operator 
for the hole with momentum k, spin a and orbital index 
a, Tik aa = d^^dkaa, and tf is the Fourier transform 
of the hopping matrix element, e a is the single-electron 
energies, and /x is the chemical potential. Introducing 
matrix notations, (ik) a p = tf an d {tn) a g = tf the 
hoppings matrix elements in the momentum representa- 



tion are given by: 

ty = 2t\ cos k 2 + 2<2 cos k 3 + 2^3 cos k\ 

+ 2i4 cos(fei + k 3 ) + 2t$ cos(k 2 + fci) + 2*6 cos(fci — k 2 ) 
+ 2i 7 cos 2fc 2 + 2t 8 cos 2k 3 + 2t g cos 2fci , (2) 

where fci = ^-k x - \k y , k 2 = k y , k 3 = %fk x + \k y . 

Note, the parameters do not change significantly upon 
changing the doping concentration. In Fig.[3]wc show two 
results of the rigid-band approximation with the Hamil- 
tonian ((T|) and the hopping values obtained in LDA calcu- 
lation for two different doping concentrations, x = 0.33 
and x = 0.61 (see Table |3|. The doping concentration 
used to calculate the chemical potential fj. was fixed to 
be x = 0.61 for both Hamiltonians. Although one finds 
the pronounced differences in the dispersion around the 
M-point, they are small around the FS. Since most of 
the physical quantities are determined by the states ly- 
ing close to the Fermi level, we can safely ignore the small 
differences of the band structure and describe the doping 
evolution of the Na 2: Co02 by simply varying the chemical 
potential. In the following we will use ab initio param- 
eters calculated for x = 0.33 and change the chemical 
potential to achieve different doping concentrations. 

Within the rigid band approximation the e' g hole pock- 
ets are well below the Fermi level for x > 0.41. Most im- 
portant however, we find the local minimum of the band 
dispersion around the T point (see Fig. [3]) to yield an 
inner FS contour centered around this point. The area 
of this electron FS pocket increases upon increasing the 
doping concentration x. As we will show later, the main 
reason for the local minimum around the T point is the 
presence of the next-nearest-neighbor hopping integrals 
which enter our tight-binding dispersion. Although this 
minimum is not yet directly observed by ARPES experi- 
ments, note that the inner FS contour would reduce the 
total FS volume and therefore may explain why the vol- 
ume of the FS observed in ARPES so far is larger than 
it follows from Luttinger's theorem 2 ^. Furthermore, an 
emergence of this pocket would influence the Hall conduc- 
tivity at high doping concentrations which is interesting 
to check experimentally. 

Note, the appearance of the inner contour of the FS 
around the T point for large doping concentrations is not 
unique to our calculations, previously it has been ob- 
tained within the LDA calculations for a single Co-layer 
per unit celU^. 



III. MAGNETIC SUSCEPTIBILITY 

To analyze the possibility of the itinerant magnetism 
we calculate the magnetic susceptibility %o(q, w = 0) 
based on the Hamiltonian Hq. The doping-dependent 
evolution of the peaks in Rexo (l, 0) is shown in Fig. 2] 
At x = 0.41 the e' g bands are below the Fermi level, and 
the FS has the form of the rounded hexagon. It results in 
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TABLE I: Single-electron energies e a (relative to e aig ) and in-plane hopping integrals t„ for Na^CoC^, where x = 0.33, 0.61, 0.7. 
(all values are in eV) 
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FIG. 3: (color online) Calculated tight-binding energy dis- 
persion, the density of states (DOS), and the Fermi surface 
for Nao.6iCo02 within the rigid-band approximation with ab 
initio parameters for x — 0.61 (the solid blue curve) and for 
x = 0.33 (red dashed curve). The horizontal (green) line de- 
notes the chemical potential li for x = 0.61. 



local minimum around the T point for the formation of 
the scattering at small momenta was originally found in 
Ref^i. 

For large x the area of the inner FS contour increases 
leading to a further decrease of the Qi. Observe that 
for x « 0.88, the FS topology again changes yielding six 
distant FS contours that moves Qi further to zero mo- 
menta. The scattering at small momenta seen in the bare 
magnetic susceptibility for x > x m is qualitatively con- 
sistent with the ferromagnetic ordering at Qfm = (0, 0), 
observed in the neutron scattering experiment a 18 i 19 ' 20 . 



IV. EFFECTS OF STRONG ELECTRON 
CORRELATIONS 



a number of nesting wave vectors around the antiferro- 
magnetic wave vector Qafm- The corresponding broad 
peaks in the Rcxo(q, 0) appear around Qafm, indicat- 
ing the tendency of the electronic system towards an 120° 
AFM SDW ordered stated. Upon increasing doping, the 
Fermi level crosses the local minimum at the T point, re- 
sulting in an almost circle inner FS contour. As soon 
as this change of the FS topology occurs, the scatter- 
ing at the momentum Qafm is strongly suppressed at 
x m > 0.56. Simultaneously, a new scattering vector, Qi, 
at small momenta appears. Correspondingly, the mag- 
netic susceptibility peaks at small momenta, indicating 
the tendency of the magnetic system towards itinerant 
SDW order with small momenta. The relevance of the 



It is important to understand the impact of electronic 
correlations on the magnetic instabilities obtained within 
the rigid band approximation. Since obtained magnetic 
susceptibility depends mostly on the topology of the FS 
one expects that the behavior shown in Fig. 2] will be 
valid even if one consider an RPA susceptibility with an 
interaction term Hi nt taken into account, at least in the 
case if the only interaction is the on-site Hubbard repul- 
sion U. The only difference would be a shift of the criti- 
cal concentrations x m , at which the FS topology changes 
and tendency to the AFM order changes towards the ten- 
dency to the FM ordered state. Similar to Refs . 16 i 30 we 
add the on-site Coulomb interaction terms to Eq. (JTJ) . At 
present, it is not completely clear to which extent the 
electronic correlations governs the low-energy properties 
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FIG. 4: (color online) The contour plot of the real part of 
the magnetic susceptibility Rexo(k,aj = 0) as a function of 
the momentum in units of 2-w/a (left), and the Fermi surface 
for corresponding doping concentration x (right). The arrows 
indicate the scattering wave vectors Qi as described in the 
text. 



in Na :c Co02 due to multi-orbital effects in this compound 
which complicates the situation. Therefore, in the follow- 
ing we discuss three different approximations valid for 
different U/t ratio. 



A. Hubbard-I approximation 

To analyze the regime of strong electron correlations 
we project the doubly occupied states out and formulate 
an effective model equivalent to the Hubbard model with 
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FIG. 5: (color online) A schematic picture of the local atomic 
states on Co and the single-particle excitations in Na^OoO^. 
Here nj, stands for number of holes, m,i enumerates single- 
particle excitations. The filling factor of the corresponding 
state upon changing the doping concentration x is given in 
square brackets. 



an infinite value of U . This approximation could be justi- 
fied by the large ratio of the on-site Coulomb interaction 
on the C0O2 cluster U with respect to the bandwidth W . 
In the atomic limit the local low-energy states on the Co 
sites are the vacuum state |0) and the single-occupied 
hole states \aa), |ei<r), \e2cr). The single-particle hole 
excitations and local atomic states are shown in Fig. [5] 
The simplest way to describe the quasiparticle excitations 
between these states is to use the projective Hubbard X- 
opcrators that take the no-double occupancy constraint 



into account automatically^: X 



\p) (q\ 



where index m «-> (p, q) enumerates quasiparticlcs. There 
is a simple correspondence between the fermionic-likc X- 
operators and single-electron creation-annihilation oper- 
ators: dfacr = J^laai^XJ 1 , where 7q CT (to) determines 

m 

the partial weight of a quasiparticle m with spin a and 
orbital index a. In these notations the Hamiltonian of 
the Hubbard model in the limit U — * 00 has the form: 

H = -J2(e p - M) Xf* -EE f fT' X f X f- ( 3 ) 

f,P f^gm,m' 

To study a quasiparticle energy spectrum of the 
system and its thermodynamics we use the Fourier 
transform of the two-time retarded Green func- 
tion in the frequency representation, G aa (k, E) = 

dj )> . This can be rewritten as: 
G aa (k,E) = E 7aa(m) 7 ^(m')^ mm '(k,S), where 

rn,m / 

t ' 



D 



(k,E) = X 



XI 



is the matrix Green 



function in the X-operators representation. 

Using the diagram technique for Hubbard X- 
operator o 32 ' 33 one obtains the generalized Dyson 
equation 3 ^: 



D(k, E) = [G^{E) - P(k, E% + S(k, E) P(k, E). 

(4) 

Here, G 1 (E) stands for the (exact) local Green func- 
tion, G% m '{E) = 5 mm >/ [E - (e p -£,)]. In the Hubbard- 
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Hubbard-I approximations one finds the narrowing of the 
bands with lowering the doping concentration x due to 
doping dependence of the quasiparticle ^s spectral weight 
introduced by the strength operator P. However, the 
doping evolution of the FS is qualitatively similar to that 
in the rigid-band picture. Namely, the bandwidth re- 
duction and the spectral weight renormalization do not 
change the topology of the FS. As a result, the presence 
of the strong electronic correlations within Hubbard-I ap- 
proximation do not change qualitatively our results for 
the bare susceptibility. Quantitatively, the critical con- 
centration x m shifts towards higher values of the doping 
and becomes x m f=a 0.68. The reason for this shift is the 
band narrowing and the renormalization of the quasipar- 
ticle's spectral weight, which enters the equation that 
determines the position of the chemical potential y,. 

Luttinger's theorem, which holds for a perturbative 
expansion of the Green's function in terms of the inter- 
action strength is violated within the Hubbard-I approx- 
imation. This violation is due to the renormalization of 
the spectral weight of the Green function by the occu- 
pation factors in the strength operator in Eq. ((4]). This 
is the reason why in spite of the e' g band narrowing the 
e' g hole pockets at the Fermi surface are still present at 
x = 0.33. 



B. Gutzwiller approximation 



FIG. 6: (color online) Calculated band structure and 
the Fermi surface topology for Na^CoCb for x — 
0.33, 0.47, 0.58, 0.68. The dashed (red), solid (blue) and dash- 
dotted (cyan) curves represent the results of the rigid-band, 
the Gutzwiller, and the Hubbard-I approximations, respec- 
tively. The horizontal (green) line denotes the position of the 
chemical potential jj,. 



I approximation the self-energy E(k, E) is equal to zero 
and the strength operator P(k, E) is replaced by the 
sum of the occupation factors, P mm (k, E) — > p mm = 

5mm' [{X p f) + (Xf 9 )] , m = m{p,q). Here '(...)' 
stands for the usual thermodynamic average. Thus, one 
obtains: 



£>(°\k,E)= Ic^-Pk 



p. 



(5) 



The Gutzwiller approximatio n 3 ^ 36 ' 37 for the Hubbard 
model provides a good description for the correlated 
metallic system. Its multiband generalization was for- 
mulated in Refi^. In this approach, the Hamiltonian 
describing the interacting system far from the metal- 
insulator transition for U >> W 



H = H 



f,a 



(6) 



with Hq being the free electron Hamiltonian {T]), is re- 
placed by the effective non-interacting Hamiltonian: 



In the paramagnetic phase the occupation factors are: 



X 



o.o 



//-•*> \ A f I ~ ~> vv 

which yields the diagonal form of the strength operator, 
P = diag x,x). Therefore, the quasiparticle bands 
formed by the a\ g — ► a\ g hoppings will be renormalized 
by the (1 + x)/2 factor, while the quasiparticle bands 
formed by the e' g hopping elements will be renormalized 
by x. 

In Fig. [5] the quasiparticle spectrum, the DOS, and 
the FS are displayed in different approximations. Within 



l-x 
2 ' 



X 







Here, lf g 

(lot, 




is the renormalized hopping, 



\-n aa i -aa — \ -0| n}ao |*0, 

bital's filling factors, x = 1 — n aa 



l faur/ 



is the or- 



is the equation for 



the chemical potential. Se aa are the Lagrange multipli- 
ers yielding the correlation induced shifts of the single- 
electron energies. The constant C is determined from 
the condition that the ground state energy is the same 
for both Hamiltonians 



(* |ffe//|*0) 



(* ff |ff|* ff > 



(8) 
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where \^q) is the wave function of the free electron sys- 
tem (J7J), and \ty g ) is the Gutzwiller wave function for the 
Hamiltonian ((6]). 

The Lagrange multipliers are determined by minimiz- 
ing the energy, 



0\H e ff |*0 



E E* 



/a d U^) +^( 9 ) 



with respect to the orbital filling factors n aa . Here C = 
^2 fe a<T n Q(T , as determined from Eq. ©. This results 

QL.a 

in the following expression for the single-electron energy 
renormalizations: 



6s c 



(10) 



It is this energy shift that forces the e' g FS hole pockets 
to sink below the Fermi energy^, which is clearly seen in 
the doping-dependent evolution of the quasiparticle dis- 
persion and the FS as obtained within Gutzwiller approx- 
imation (Fig. [5]) • Although the narrowing of the bands 
due to strong correlations is similar to the one found in 
the Hubbard-I approximation, the FS obeys Luttinger's 
theorem. Note, in contrast to the Hubbard-I approxi- 
mation the relative positions of the <2g-bands are also 
renormalized by 8e aa . 

At the same time, for x > 0.4 the topology of the FS in 
the Gutzwiller approximation is qualitatively the same as 
in the rigid-band picture. The also yields similar results 
for the bare susceptibility's doping dependence discussed 
in Section IIIII The only effect of the strong correlations 
for xo is the observed shift of the critical concentration 
towards higher values, x m 0.6. This is due to combined 
effect of the bands narrowing and the doping dependence 
of the aig and e' g band's relative positions, determined by 
the Eq. ([TUj). 

Note, for x < 0.4, due to different FS topology that 
occurs in the Gutzwiller approximation, the bare sus- 
ceptibility differs from that obtained in Refi^ where the 
strong renormalization of the electronic bands removing 
e' g pockets away from the FS was neglected. 

C. FLEX approximation 

A certain disadvantage of the Gutzwiller and Hubbard- 
I like approximations is that the dynamic character of 
electronic correlations is not taken into account within 
these approaches. At the same time, the momentum and 
frequency dependencies of the self-energy S(k, u) play 
a crucial role, in particular, for determining the low- 
energies excitations close to the Fermi level. In this sub- 
section we focus on the ai 5 -band with nearest and next- 
nearest hopping integrals only and employ the single- 
band Fluctuation Exchange approximation (FLEX)^ 



which sums all particle-hole(particle) ladder graphs for 
the generating functional self-consistently valid for the 
intermediate strength of the correlations. The FLEX 
equations for the single-particle Green function G, the 
self-energy E, the effective interaction V, the bare (x°) 
and renormalized spin (x s ) and charge (x°) susceptibili- 
ties read 



G k (w n ) 



[W„ - Tk + jJL - Ek(0>n)] 1 



mi 



T 



Jr ^ Vk- P (w„ - w m )G p (w m ), (12) 



p,m 



Vq(fm) = U 2 



^Xq(^m) + ^Xq(^m) - XqOm) 



(13) 



X>m) 



Xq C ("»»») 



T 

'N 



^ G k+q (w„ + C m )GkK), ( 14 ) 



x> m ) 



iTUxlM 



(15) 



where Lu n = i j nT{2n + 1) and v m = ii:T{2m). In the last 
equation the '— ' sign in the denominator corresponds to 
the Xq(^m), while the '+' sign corresponds to the XqC^m)- 
We compute the Matsubara summations using 'almost 
real contour' technique of Ref4£. I.e., the contour inte- 
grals are performed with a finite shift (0 < 7 < iTn/2) 
into the upper half-plane. All final results are analyti- 
cally continued from w + 27 onto the real axis to + i0 + by 
Pade approximation. The following results are based on 
FLEX solutions using a lattice of 64x64 sites with 4096 
equidistant w-points in the energy range of [—30, 30] . The 
temperature has been kept at T = 0.05r, where r is the 
hopping amplitude to the nearest neighbors for the a\ g 
band corresponding to t" ls > ai! > = t^ 13 = t" 1 ^ 1 '. The 
Hubbard repulsion was set to U = 8r. 

Previously, the FLEX approximation has been applied 
successfully to the study of superconductivity as well as 
spin and charge excitations in Na^CoO i 41 ' 42 . Comple- 
mentary, we will focus on the quasiparticle dispersion and 
study the impact of the momentum and frequency depen- 
dencies of the S(k, ui), and the role played by the next- 
nearest hopping integral, r', corresponding to t^ lg 19 = 
^ai g a lg _ j_a lg a lg ^ quasiparticle dispersion E^, which 
is determined from equation — Tk + (i — Sk(^k) = 0, 
is shown in Fig. [7] for r' = and r' = —0.45, in units of 
t. First, observe that the local minimum around the T 
point appears only if the next-nearest-neighbor hopping 
t' is included which agrees with our previous findings. In 
addition, we obtain a pronounced mass enhancement of 
the order of unity at the FS crossings - the so-called kink 
structure. This enhancement is due to low-energy spin 
fluctuations which are present in Xq( w )^- 

To shed more light onto the two-dimensional spin cor- 
relations, in Fig. [5] we display the static spin structure 
factor Rex£.(ix> = 0) from the FLEX for two different 
doping concentrations. As doping increases from x = 
towards x = 0.35 the maximum in the spin susceptibil- 
ity Xk( w = 0) moves towards the K-point of the first 
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FIG. 7: Quasiparticle dispersion (in units of r, relative 
to fi) within FLEX approximation for (a) r' = and (b) 
t' = —0.45 and for two doping concentrations. 



r km r 

FIG. 8: Doping dependence of the static spin structure factor 
Rex k (w = 0) for (a) r' = and (b) r' = -0.45. Notice that 
for large U = 8t the commensurate peak at K point is absent 
at a very low x. 



BZ and develops into a sharp and commensurate peak 
at Qafm and the incommensurate spin fluctuations are 
suppressed. One may also note that the commensurate 
peak is ~ 60% larger for r' = —0.45 than for r' = 0. 
These results are consistent with those obtained in a pre- 
vious sections. Wc further notice smooth evolution of the 
quasiparticle dynamics with doping in Na^CoC^ showing 
no sign of unusual behavior at x = 0. 

The frequency dependence of the imaginary part of 
the quasiparticle self-energy, i.e. ImEk) near the FS is 
shown in Fig. [5J We find the self-energy to be nearly 
isotropic along the FS with only a weak maximum oc- 
curring into the direction of the commensurate spin fluc- 
tuations. Near the Fermi energy the self-energy is clearly 
proportional to lo 2 at low energies for all dopings shown, 
which is indicative of the normal Fermi-liquid behavior. 
This is in sharp contrast to the FLEX analysis of the 
Hubbard model on the square lattice close to half-filling. 
There one typically finds 'marginal' Fermi-liquid behav- 
ior with ImEk ~ u> over a wide range of frequencie s 43 ' 44 . 
Therefore, along this line one is tempted to conclude that 
the normal state of the superconducting cobaltates is 
more of conventional metallic nature than in the High-T c 
cupratcs. This is even more so, if one realizes from Fig. [9] 
that the quasiparticle scattering rate displays its smallest 
curvature for x = 0.35, which implies the quasiparticles 
to be rather well defined there. For lower x proximity 
of the FS to the van Hove singularity (see flat region of 
dispersion in Fig. [7]) enhances both the absolute value of 
ImEk oc a; as well as the curvature. This effect is most 
pronounced for r' = 0. 
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X=0.35 





2 1 -1 CO 10—10) 



FIG. 9: Frequency dependence of quasiparticle self-energy 
—ImEk near the FS in direction T — K for (a) r' = and 
(b) t' = -0.45. 



V. CONCLUSION 

To conclude, we have calculated the doping dependent 
magnetic susceptibility in the tight-binding model with 
ab-initio calculated parameters. We find, that at a crit- 
ical doping concentration, x m , electron pocket develops 
on the FS in the center of the Brillouin zone. For x < x m , 
the system shows a tendency towards an 120° AFM or- 
dered state, while for x > x m a peak in the magnetic 
susceptibility forms at small wave vectors indicating a 
strong tendency towards an itinerant FS state. Within a 
tight-binding model we have estimated x m to be approx- 
imately 0.56. Analyzing the influence of strong Coulomb 
repulsion and the corresponding reduction of the band- 
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width and the quasiparticle spectral weight in the strong- 
coupling Hubbard-I and Gutzwiller approximations, we 
have shown that the critical concentration changes to 
x m w 0.68 and x m ~ 0.6, respectively. At the same time, 
the underlying physics of the formation of the itinerant 
FM state remains the same. 

We neglected the bonding-antibonding splitting due to 
the 3-dimensionality in the non- intercalated compounds. 
This splitting was taken into account in Ref4^, where 
within the FLEX approximation the single a lff -band 
Hubbard model was considered. The results obtained 
also suggest a tendency to FM fluctuations for high dop- 
ing concentrations. The presence of a local band mini- 
mum around the T point played a crucial role, similar to 
our present study. 

To analyze the low-energy quasiparticle properties at 
low doping concentrations we have employed the single- 
band Hubbard model within the FLEX approximation. 
We have found a significant FS mass enhancement of or- 
der unity due to quasiparticle scattering from spin fluctu- 
ations. In contrast to the Hubbard model on the square 
lattice we have found the quasiparticle scattering rate to 



display a conventional Fermi-liquid type of energy depen- 
dence. We have also shown that the static spin structure 
factor exhibits a large commensurate peak at wave vec- 
tor Qafju for doping concentrations of x ~ 0.35. This 
response was found to be significantly enhanced by the 
next-nearest-neighbor hopping, emphasizing its signifi- 
cance. 
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